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ABSTRACT 


We examine the Galerkin (including single-mode and Lorenz) 
equations for convection in a sphere to determine which physi- 
cal processes are neglected when tht equations of motion are 
truncated too severely. We test -<ar conclusions by calculating 
solutions to the equations of motion for different values of 
the Rayleigh number and for different values of the limit of 
the horizontal spatial resolution. We show that the transitions 
from steady-state to periodic, then to aperiodic convection 
depend not only on the Rayleigh number but also very strongly 
on the horizontal resolution. All of our models are well- 
resolved in the vertical direction, so the transitions do not 
appear to be due to poorly resolved bo undary- layers . One of the 
effects of truncation is to enhance the high wavenumber end of 
the kinetic energy and thermal variance spectra. Our numerical 
examples indicate that as long as the kinetic energy spectrum 
decreases with wavenumber, a truncation gives a qualitatively 
correct solution. 
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In Rayleigh- Bernard convection, discrete transitions from 
steady-state to periodic to aperiodic convection have been experi- 
mentally observed. (See the recent reviews by Fenstermacher et al. 
1978 and Basse 1978.) As the Rayleigh number is increased and the 
fluid becomes more "turbulent”, the Fourier spectrum (in time) of 
the velocity develops a single spike (and its overtones) and 
shows a gradual increase of the broad band background noise 
that eventually overwhelms the spikes. Although the transitions 
depend not only on the Rayleigh number but also on the 
Prandtl number and initial conditions, there has recently been 
much interest in trying to compute these transitions 
from the actual equations of motion. 

In attempting to compute time-dependent numerical solutions 
to the three-dimensional Navier-Stokes equation, one is forced to 
make severe approximations. When simplifying the equations of 
motion to make them numerically tractable, one hopes to esta- 
blish a compromise so that the modified equations are uncompli- 
cated enough to be easily solved, yet complete 
enough that the underlying physics of the fluid dynamics is 
not lost. 

The crudest approximation is the Lorenz (1963) model. The 
Lorenz model predicts not only the transitions to steady-state 
and time-dependent convection, but also a sequence of bifurca- 
tions that eventually leads to chaotic (aperiodic) behavior. 

For low Rayleigh numbers near the onset of convection the heat 
flux (Nusselt number) predicted by the Lorenz model is in fair 


agreement with laboratory results. As the Rayleigh number increases, 
the calculated and experimentally observed Nusselt numbers begin 
to differ. When the Rayleigh number is as large as the one at 
which the Lorenz model predicts a transition to chaos, there is 
an appreciable difference between theoretical and experimental 
values of the heat flux and one must seriously question the quali- 
tative behavior of the time-dependency of the solution. McLaughlin 
and Martin (1975) have expanded the Lorenz model to four inter- 
acting modes and have found support of the Rouelle-Takens (1971) 
theory of turbulence, which states that after no more than three 
bifurcations to a periodic or quasi-periodic state there should 
be a transition to aperiodicity. The fundamental question to be 
answered, of course, is whether the qualitative time dependence 
of these equations is due to the underlying physics that these 
equations are trying to model or whether bifurcations are a 
general property of sets of severely truncated nonlinear differen- 
tial equations. A truncation of the governing equations of con- 
vection that is less severe than McLaughlin and Martin's treat- 
ment in the radial direction is single-mode theory (Gough et al . , 
1975). Single-mode theory has only one horizontal mode 
so has less horizontal resolution than McLaughlin and Martin's 
4-mode solution. Surprisingly, the numerical solutions to the 
single-mode equations (Toorare et al., 1977) do not exhibit bifur- 
cations to periodic or aperiodic states and are time-independenx 
for all Rayleigh numbers. Numerical solutions to a truncated 
Galerkin expansion of the equations of convection that is less 
severe than both single-mode theory and McLaughlin and Martin's 
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equations have been computed by this author (Marcus 1980a, 1980b), 
in which not just one but several horizontal modes of the expan- 
sion are retained. For steady-state convection the results are 
in good agreement with single-mode calculations. However, at lar- 
ger Rayleigh numbers we find that the solutions become periodic 
in time and, as the Rayleigh number is increased further, aperio- 
dic in time. For some Rayleigh numbers that produce steady- 
state solutions, we find (holding the Rayleigh number and the 
resolution in the radial direction fixed) that as we decrease 
the number of horizontal modes in the Galerkin expansion, there 
is a transition from steady-state convection to a solution that 
is periodic in time. As the number of modes is decreased still 
further, the solutions become aperiodic. Obviously, the bifur- 
cations depend not only on the Rayleigh nximber, Prandtl number 
geometry and initial conditions , but also on the horizontal 
resolution of the equations of motion. 

In trying to understand mathematically the bifurcation 
sequence of a truncated representation of the equations of motion, 
it is easy to lose sight of what is physically happening in 
the fluid. Therefore, the purpose of this paper is to examine 
the solutions to truncated modal equations for convection in a 
sphere and to determine which qualitative features of the solu- 
tions represent real physical processes in the fluid and which 
features are due solely to the effects of truncation. 

In section 2 of this paper we briefly review the Galerkin 
multi-mode equations (including single-mode and Lorenz) for 


spherical convection. We attempt to describe the physics that 
each system of equations models, which physical processes are 
neglected by the various truncation schemes, and what artifi- 
cial constraints each model imposes on its solutions. In the 
third section we present the results of our multi-mode calcula- 
tions for two different values of the Rayleigh number. For 
each Rayleigh number we compute several models, each with a 
different degree of horizontal truncation. By computing how 
the energy spectra, convective flux, and temperature gradient 
change as a function of the severity of truncation, we provide 
a possible explanation for the time-dependence of our solutions. 
Our conclusions appear in section 4. 

II. APPROXIMATIONS NEEDED FOR THE LORENZ, 

SINGLE-MODE AND MULTI-MODE MODELS 

Convection in a Boussinesq fluid is governed by the Navier- 
Stokes, continuity and thermal diffusion equations, and the 
Boussinesq equation of state. (See, for example, Chandrasekhar , 
1961). A standard technique used to simplify these coupled, 
nonlinear, partial differential equations is the Galerkin method. 
The thermodynamic quantities and velocity are expanded as an 
infinite sum of coefficients multiplied by orthonormal functions 
and substituted into the governing equations. Then, depending 
on how many of the coefficients are solved and how many are 
arbitrarily set equal to zero, one arrives at a Lorenz, single- 
mode, or m.ulti-mode model. 


i 
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A) Review of the Multi»Mode Equations 


Let us consider convection in a self- gravitating sphere of 
Boussinesq fluid with thermal expansion coefficient a, heat 
capacity Cp, kinematic viscosity v, thermal diffusivity k, 
radius d, ajid a heat source H(r) in the fluid. Each 

scalar quantity, such as the temperature, is written as a 
sum of its mean, <T(r,t)>, and fluctuating, T(r ,9 ,t) , parts where 


<T(r 


,t)> = -/* T(r ,6 ,(}) ,t) 


dfi/4TT 


( 2 . 1 ) 


and 


47T 


T (r,6,0,t) =2(27r)^/^ 


T. Y. 


m-.l 


* I" 

Re(Y^’^) and Im(Y^’”^) are the real and imaginary parts of the 
of the spherical harmonic. The velocity is written as a sum 
of its poloidal Vp and toroidal v^ parts which are derived from 
scalar fields u and ^ 

V = 7 C3(ru)/3r] - (rV^w) i (2.3) 

~p r 



v,p s r7 X (\pi^) 


(2.4) 


Substituting expressions (2.3) and (2.4) into the equa- 
tions of motion yields the equations for the coefficients for 
the temperature, T, pressure, ?, gravitational potential, 
and velocity, v (Marcus, 1980a): 


/ at 


Su 


Y , i,ns 


-r C4CU1)]-1 [R,PrrT^_j_„ ♦ 3 


+ ?r ©j '“Y.t.m’ 

- lUA+l)"^ {r • C(v7)v]} . 


m 


-2(U1)-^ (r • 7 X Cv-7)v:)^_j_^^ 

3 T„ « . / 37 * 1 


(2.5) 


( 2 . 6 ) 


- aU+l) 0<T>/3r) w . _/r 

Y , z ,m 


- [v • 7T] 


Y,z,m 


- r"^3 {rCr 

. {7 . C(v . V)v]}^^,^^ 
'* Y * Z , n ^ " 


(2.7) 


( 2 . 8 ) 

(2.9) 


< y > S<'^>sO 

s { 3 (r" 3< T>/ Sr) / 3r + S*^/3r 

-3 [/ " i (£+ l ) T , u 0 _] / 3 r } 

-d L Z-. - ■^'‘ y,i^m Y>Z,m 

Y ) Z ) n, 


(2 .10 ) 


( 2.11 ) 
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where is the differential operator defined by its action on 


the scalar » f: 


(f) ! [9*(rf)/8r^ - HU+1) f/rV: 


i 

and where«^(r) is the luminosity 


^(r) = 4 tt J <H> r ^ dr’ 


( 2 . 12 ) 


(2.13) 


In equations (2.5) - (2.11), Rs = aGd^i,(d)/3k^vC„ is the Rayleigh 

P 

number, Pr s v/k is the Prandtl number and y stands for either 
R or I. 

2 

In equations (2.5) - (2.13) the unit of time is k/d , 
length is d, mass is pd^ and temperature is jC( 

Equations (2.5) - (2.10) may be thought of as the governing 
equations for each eddy or mode (Y,i,m) that make up the total 
velocity field. The nonlinear terms in equations (2.5) - (2.13) 
such as (re^'Cv* V)v]> . are the eddy-eddy interaction terms, 

with contributions from all other pairs of modes 
and (y" , m” ) that obey certain selection rules. The selection 
rules and the explicit expressions for the nonlinear inter- 
actions are given in a previous paper (Marcus 1979) in terms of 
Wigner-3j symbols. For a sphere with an impermeable, stress- 
free boundary the velocity is constrained at r=l so that: 


w c - (1) = 0 

f j 
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(2. lii) 
(2. 15) 


0 


(2. 15) 


lo- 


ws Also rsquirs that the surface be isothermal: 


^Y,i,m ° • 


(2.17) 


We are free to choose the mean temperature to be zero 

at rsl: 


<T(r»l)> s 0 . (2.18) 

However, the gradient of the mean temperature (and therefore 
the flux) at rsl is free to vary. The central temperature, 

<T(0)>, is also free to vary and is a measure of the efficiency 
of the overall convective flux. The lower the value of <T(0)>, 
the more isothermal the fluid is. The central temperature is given 


<T(0) 


> = - feC./r^ dr 

I 

- I V JUi,+l) T 
j Y.i.m 


, w . /r dr 

Y,i,m 


f [f r'^ <T(r’)> dr’]/r^ dr 
*0 -^0 


(2.19) 


We must use the central temperature as a measure of the effi- 
ciency convection because the Nusselt number is not well-defined 
for our boundary conditions. 


3) Sufficient Conditions for a Good Truncation 


The inrinite se 


t modal ecuaticns (2.5) - (2.13) for the 


coefficients can only be solved by arbitrarily setting some of 


■one coemcients ecual 


to cero (or some other f'cnctional form) 


and extlicitly solving fcr the re.maining finite set of cceffi- 


REPRODUciBn.rry op the 

ORIGl^JAT. PAGE IS POOR 

cients. What are the consequences of setting some modes equal 
to zero? Thte equation for mean value of the temperature, 
(2.11), is well approximated, if and only if the term^iU+l) 
“Y»i,m summed over the finite set of )<ept modes, 

is nearly equal to what it would be if it were summed over all 
modes. Now, ViU+l) w . T * /r is equal to the con- 

vective flux and the contribution from each mode is just the 
convective flux carried by that particular eddy. Therefore 
equation (2.11) is well-approximated if we keep those 

eddies that carry most of the flux in the Galerkin expansion. 
Similarly it can be shown that equations (2.5) - (2.8) are well 
approximated only if we include the modes that are responsible 
for (1) the production of kinetic energy from buoyancy forces, 
(2) the production of the temperature variance, 1/2 T , (3) the 
viscous dissipation of kinetic energy, ( 4 ) the dissipation of 
the temperature variance, (5) and those modes that provide 
the nonlinear cascade of energy from the production modes to 
the dissipative modes. We expect that the modes most responsible 
for production of the kinetic energy temperature variance and 
convective flux are the largest spatial modes. We also expect 
that if we wish to include all of the modes that are important 
in the cascade and dissipation of kinetic energy and temperature 
variance, we will have to retain all modes with Reynolds or 
-eclet numbers are greater than one. 
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C. Tha Effects of Truncations on the Kinetic Energy 

The rate at which kinetic energy-y|— J y v^d^r enters 
the fluids due to buoyancy is (Marcus, 1980a) 




47t Pr Rs P ^ 

J fl 


Y.A.Jn 


„ Wy 1 , _ r" dr (2.20 
Y » * ,n Y ** *® 


There are no cross-terms between different modes on the right- 
hand side of equation (2.20) and each term represents the kinetic 
energy contribution from one mode (t,m,a). However, combining 
equation (2.11) with (2.20) shows us that we can write in 

terms of the luminosity and temperature gradient: 


KE. = 4ff Pr Rs 
in 


r 

0 L 


5<T> 


.r=Z! 


r 

- ^ r J r'^ <T(r')> dr’ 


dr 


(2.21) 


By numerical experimentation we have found that no matter 
how few modes are kept in the Salerkin expansion, the mean temp- 
erature gradient becomes nearly isothermal in the sense that 




(2 


) 


Using equation (2.22) and taking the time-average (denoted by 
double brackets) of equation (2.21) we obtain: 


<<KE. >> 
in 


4ir ?r Rs /X (r) rdr 
0 


( 2 . 2 ?) 


We find that even the most severe truncations produce a close 
approximation to the correct value of <<?(E. >>. 
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The time-averaged value of the rate at which kinetic 
energy is dissipated » must be equal to <<KE^^>>. 

is given by 


KE^Ut = 


y >-Ita+l)3^ r"^ 


Jt,m 


-2 


-JUZ+1) r ‘ 3(ru> 


J^H r"^a^(rKE)/3r^ (2.24) 

)/3r 3Cr J). (w 


Y,t,m Y»^»m 


)]/3r 


2 

r^dr 


where KE(r) is the kinetic energy of the fluid at radius r 
and is 


KE(r) 


I 


UJl+l) 


y,i,m 






(2.25) 


Again, there are no cross terms between modes on the right- 
hand side of equation (2,24) and each term in the sum repre- 
sents the dissipation due to one mode. If the high wavenumber 
modes responsible for the viscous dissipation are not inoluded 
in the Galerkin expansion (or if the modes that are responsible 
for the cascade of kinetic energy to the dissipative modes are 

not included) <<KE. >> will net be stronglv affected. However. 

in 

to keep <<KE >> equal to <<KE- >> the fluid must compensate 
^ out ^ in 

by dissipating more kinetic energy in the large scale modes. 
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Trom equation (2.24) we see that one way in which the rate of 
dissipation can be increased is by increasing the kinetic energy 
of the modes. We therefore expect the kinetic energy of a se- 
verely truncated system to be abnormally high. This increase 
will be evident in the numerical examples in the next section. 


D. The Effect of Truncation on the Fluctuating Thermal Energy 


The rate at which temperature variance is created in the 
fluid is 


j 

TE. = -47t / 

Jn 


3<T> 

9r 


I 

Y,^,m 


T , . iL(il+l) dr (2.27) 

y,i,m Y,£,m 


Each term in equation (2.27) corresponds to the thermal 
input of one mode. Even though | ^ ■ [ will generally be much 
less than w® have found that ror fixed Prandtl and 

Rayleigh numbers, I'^^l vary by an order of magnitude 
depending upon the number of modes kept in the Galerkin expan- 
tion. Therefore, <<TE^^>> (unlike <<KE^^>>)is a sensitive func- 
tion of the truncation. The rate at which the temperature vari 
ance is dissipated is 


I 


TE^ ^ = -4tt L. + iUil+l) r"^T^ . ^3 r^dr 

out Y»^»m Y>Jt,m 

y 5 X 


(2.28) 


If the Galerkin truncation does not include the thermally . 

dissipative modes, the truncated solution will have to adjust 

itself so that <<TE .>> is kept equal to <<TE. >>. The solu 

out in 
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tion can increase the rate of thermal dissipation in the re- 
tained modes by increasing the thermal variance of the modes. 

However, unlike <<KE. >>, <<TE,. >> is not constrained and the 

in in 

fluid can adjust to its inability to dissipate the thermal 

variance by decreasing <<TE. >>. Since <<TE. >> is propor- 

in in 

tional to the mean- temperature gradient (eq. 2.27), the fluid 
can reduce its rate of production of thermal variance by be- 
coming more isothermal. In the next section we show numerical 
examples in which a truncated solution both increases 
by increasing its thermal variance and decreases by 

becoming more isothermal. 

E. Single-Mode Theory 

The severest truncation of a multi-mode expansion is to 
retain only one horizontal mode. This requires that the solu- 
tion be of the form: 


T(r,9,4»,t) = <T(r,t)> + T(r,t) h (6,4)) (2.29) 

P(r,9,4),t) = <P(r,t)> + P(r,t) h (6,4') (2.30) 

oj(r,6,t) = ai(r,t) h (0,4>) (2.31) 

4; = 0 (2.32) 


where h(e,4)) is an eigenfr.nction of the horizontal Laplacian, 

(v' - 

Because the toroidal modes are not involved in the convective 
flux, kinetic energy production, cr temperature variance produc- 
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tion they are neglected in single-mode theory. Our multi-mode 
numerical e:cperiments have shown that the toroidal velocity is 
much smaller than the poloidal velocity except for large wave- 
number modes in large Rayleigh number convection. (See 
Marcus 1980b). 

Unlike expansions with more than one horizontal mode, the 
single-mode solutions are always time-independent. Toomre et al. 
(19 77), working with a plane-parallel geometry> also found that 
a single-mode always leads to a steady-state solution. Expansions 
with a single-mode suffer not only from the effects of truncation 
mentioned in the previous section, but also from other problems. 
For example, the correlation between the radial velocity and 
temperature , 


'V ^ 1 7 1/9 

6 = <TV >/<T >“ <v > 
r r 


(2.33) 


is always identically equal to 1 for a single-mode; whereas, 

experimentally, Deardorff and Willis (1967) have found that the 

correlation in air for Rayleigh-Bernard convection is between 

5 7 

.5 and .7 for Rayleigh numbers between 6><10 and 10 . The con- 
vective flux, <TV^> , that is far from the bounday predicted by 
single-mode theory is in good agreement with the flux predicted , 

multi-mode calculations (see §3). Because the single-mode over- 

'^2 

estimates 6, it always underesti.mates <T ><V^>, the product or 
the thermal variance and radial component of the kinetic energy. 
Another peculiarity of the single-mode equations is that the 
thickness of the boundary -layer at the surface i' conirolled by 
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viscosity and decreases as the Rayleigh number is increased 
(see Toomre et al. 1972) In a real fluid we would expect the 
boundary-layer to become turbulent and wide as the radially 
moving fluid smashes into the impermeable outer boundary. The 
thickness of the turbulent boundary-layer is not regulated by 
viscosity, but by the rate at which energy can be transferred 
to other modes. The increase in boundary-layer thickness due 
to the nonlinear cascade in a multi-mode calculation has been 
reported by this author elsewhere (1980a). In a single-mode 
calculation with a large Rayleigh number and an artificially 
thin boundary-layer, most of the dissipation of kinetic energy 
takes place near the surface with 

<<K£ .>> ~ 4ttP / (v • 7^ v) r^ dr (2.34) 

out r / — — 

•'i-x 

where X is the thickness of the boundary-layer. From equation 
(2.34) we see that proportional to 1/X. Therefore, 

a single-mode calculation can compensate for its loss of dissi- 
pation in the missing high wavenumber modes by decreasing X. 


F. Lorenz Model 

A further truncation of single-mode expansion gives us the 
Lorenz model. Using the equilibrium conductive temperature 
gradient with the single-mode equations, we can compute the com- 
plete set of orthonormal eigenmodes of the velocity and 
temperature (as functions of radius). By expanding the radial , 
dependence of the velocity and temperature in terms of these 
eigenmodes, substituting the expansions into the single-mode 
equations and retaining only a single mode in the radial expan- 
sion, we obtain the Lorenz equations. These equations 
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wsr« originally derived for a convecting fluid in a plane-parallel 
geometry, but they can easily be extended to a spherical geometry. 

The Lorenz model not only suffers from all of the physical approxi- 
mations of single-mode theory but also contains some additional lia- 
bilities. Because the functional form of the velocity and fluc- 
tuating temperature are fixed and only their amplitudes are allowed to 
vary, the fluid can never develop boundary- layers to help dissipate 
the kinetic and thermal energy. More importantly, becauta the 
functional form of the velocity and temperature are fixed, the 
mean-temperature gradient can not become isothermal. 

If we were interested in computing solutions only when the 
Rayleigh number is slightly greater than its critical value, it 

would be practical to expand the velocity and temperature in the 
eigenmodes that are calculated with the conductive temperature gradient . 
However, these are .not a very useful set o.f functions in which to 
expand the velocity and temperature when the Rayleigh number is 
large. For example, for any large Rayleigh number, we can choose a 
complete basis in which to expand the velocity and temperature by 
calculating the fundamental and all of the higher harmonic solutions 
to the single-mode equations. By retaining only the fundamental 
mode in the expansion, a modified set of Lorenz-type equations is 
obtained. We have computed the steady-state solutions to the 
regular spherical Lorenz equation and to the new modified Lorenz 
equations for a Rayleigh number 30 ti.mes greater than “he criti- 
cal value for the onset of convection. The solution to the regular 


f Lorenz equation is unstable with respect 

i 

r . . 

i bations ; both the solution to the mcdcrc 


to ti.me-dependent pertur- 
ed Lorenz equation and 
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the steady-state solution to the multi-mode equations with this 
Rayleigh number are stable. We conclude that qualitative descrip- 
tion of the Lorenz model is not accurate for large Rayleigh numbers. 

III. NUMERICAL RESULTS OF MULTI-MODE CALCULATIONS 

In this section we present the numerical results of 

multi-mode calculations for Pr=10 and Rayleigh numbers of 
4 5 

10 and 10 . For each Rayleigh number we repeat the cal- 
culation several times, each time using a different set of 
modes to show the effects of truncation. For all calculations, 
the heat source H(r) (see eq.2.13) is constant for riO.3 and 
zero elsewhere. 

A. Rs = 10^^ 


’ To compute solutions to the modal equations, we have chosen 

1 ^^ , the set of modes in the Galerkin expansion to be all of the 

spherical harmonics, with I i The radial 

dependence is finite— differenced with 128 grid points. For 
^cutof f ~ ^ ^ ^ find that the solution is time-indepen- 
dent. A complete description of the solution with I = 12 

cutoff 

appears elsewhere (Marcus 1980a). To compare the overall fea- 
tures of the truncated solutions, we have listed the central 

temperature, KE . and TE. as a function of I * in Table 
XU xn cutorr 

3.1. 


TABLE 3.1 


^cutoff 

<T(0)> 

KE. 

in 

in 

12 

0,666 

4.94x10^ 

CO 

o 

4^- 

9 

0.685 

4,94x10^ 

4.01 

5 

0.675 

4.90x10® 

c 

3.93 

3 

0.528 

U . 44x10 

1.71 
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There is virtually no difference in the calculated values of 
<T(0)>, or for ^cutoff"® 12, -which indicates 

that modes with i > 6 are not important in production, trans- 
port or dissipation of energy. Using the value of KE^j^ from 
table 3.1, we find that the Kolmogorov length is 'v 0.212 which 
approximately corresponds to a wavenumber, i, of 'v 4. The 
solution with shows the effects of truncation; the 

rate of input of thermal energy for is nearly 60% 


lower than it is for ^cutoff"*^ 


2. The rate KE. for 

in cui.o^i. 

is nearly equal to KE^^^ for ^cutoff large decrease in 

TE. is consistent with the analvsis presented in § 2 which shows that the 

fluid can compensate for the loss of the thermally diffusive 

modes bv decreasing TE. . KE-„ is constrained by the fact that 

it must always be approximately equal to 4TrPrRs / Xi’r)r (ir = 4 . 94x10^ . 


To compensate for the loss of the high wavenumber modes that 

dissioate the thermal variance when l * = 3 , the fluid de- 

cutorr 

creases TE. by making the temoerature gradient more nearlv 
isothermal. The isothermal nature of the ^^utof**^ ” ^ solution 
can be seen by noting that the central temperature for l 


cutorf 

= 3 is less than it is for ^^;^.toff ” 

A more sensitive probe of the effects of truncation is the 
kinetic and thermal energy spectra as functions of the horicon- 
tal wavenumber. In table 3.2 we have listed TE(Z ,r=0 . 5 ), which 
is the 2-dimensional thermal variance spectr’cm at r=0.5, with 
wavenumber i , i . e . 


f ,m 


(3.1) 
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We have also listed the kinetic energy spectra, K£(£,0.5), at 
r * 0.5, as functions of l and ^^y^off 3.3. The ki- 

netic energy spectra show that the value of ^^^cutoff 
is higher than it should be. As pointed out in 1 3, the trun- 
cation causes an upward curl in the energy spectrum at 
because the energy that cascades down from the large scale 
modes piles up at The upward curl at the large wave- 

number end of the spectrum is even more pronounced in the thermal 
variance spectra. Because the Prandtl number is greater than 
unity, the dissipation of thermal energy is less efficient than 
the diffusion of kinetic energy. The thermal variance does not 
dissipate in the production modes as does the kinetic energy 
and is free to cascade down the spectrum and pile up at the 
large wavenumbers. For the severest truncation, ® 3, 

the thermal energy spectrum has inverted itself and TEC 3, 0.5) > 
TE(2,0.5) > TE(1,0.5). 

B. Rs = 10^ 

For a Rayleigh number of 10^ and a Prandtl number of 10, 
we have computed solutions for ^g^^gff ® 12,9,6,4,3,2, and 1. 

With *-gu^gff ® 1 the solution is steady-state and the multi-mode 
equations reduce to those of single-mode theory. For a compari- 
son between single and multi-mode solutions, we have plotted 
the kinetic energy of the i - 1 mode as a function of radius 
in figure 1 (solid lin*.). Superimposed on this figure is the 
kinetic energy of the £ = 1 mode (broken line) computed from the 
steady state solution of the multi-mode equation with . 


12. The functional form of the two curves is quite similar, 
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th« main dlfftrsnct being thet the eingle-mode klnetie ener^ 
ie consietently higher then the multi-mode solution. This 
difference in height confirms the predictions we made in 12: 
the kinetic energy of the single-mode must be enhanced to in- 
crease its rate of viscous dissipation. For the single-mode* 
KEout ^ whereas for the l • I component of the 

g 

multi-mode solution,KE^^^ is only 3.13 x 10 . Approximately 
321 of the kinetic energy produced in the i « 1 component of 
the multi-mode solution is lost not through dissipation but 
through the nonlinear energy cascade. 

In figure 2 we have plotted the temperature variance of 
the i « 1 mode of the multi-mode solution (broken line) and 
the single-mode solution (solid line). As in figure 1, the 
two curves have the same function form* but, in general, the 
single-mode thermal variance is greater than the multi-mode 
variance. The greater thermal variance allows the single-mode 
to increase its rate of thermal dissipation. The rates at 
which the temperature variance is dissipated from the £ s i 
component of the single- and multi-mode solution are 0.293 and 
0.231 respectively. 

For all solutions computed with ^^y^off - solutions 

are steady-state and show truncation effects similar to those 

4 

found for Rj s 10 . For ~ *♦* the temperature spectrum is 

inverted with TE(l+l, 0.5) > TE(l, 0.5). The kinetic energy 
spectrum is not inverted. In figure 3 we have plotted 


Q i K£ Ub2, rsO.S}/ KE(t«3, r«0.3) at a ftir.w'rion 
of Q is a maasura of tha upward curl of tha kinatic 

anargy spactrua a- ta3. If thara wara no truncation af facts, 
wa would axpact Q always to ba graatar than 1. Zf Q bacomas 
lass than l,it aaans that tha kinatic anargy spactrua is 

invartad. i.a., X£U>3, r«0.5) > K£Cts2, rsO.S). Figura 3 
shows that Q is graatar than 1 but dacraasas as i^utoff ^*^^*^***‘ 
By axtrapolating tha points in figura 3, wa may axpact that Q :.s 
lass than 1 for ^gytoff*^* ^cutoff*^ solution is no 

longar staady-stata but is periodic in tiaa. Tha kinatic anargy 
calculated with * 3 at r*0,5 as a function of wavelength, 

i, and as a function of tiaa is plotted in figur-a u for one 
period of tha fluid's oscillation. 

Wa have ar' . trarily labelled tha left-hand axis of f igura 
^ as t»0 but,infact, it takes many iterations for the transients 
in the fluid to settle down and for the motioris to become per- 
iodic. At tsO, the kinetic energy of the 1*1. 2 and 3 wavelengths 
are similar in value to the stationary values obtained with 
tcutoff*^^* increases, tha kinati.? energy of 1*2 and i*3 

modes increases ; they are unable to dissipate th</;ir kinetic anargy 
as fast as it cascades into (or is produced in) 
the modes. At t*.0467 the kinatic energy of the J*1 modes 
becomes lass than that of the i*2 mode, and at 
t*.0503 the kinetic energy of the t*l and 1=3 modes cross. At 
this point in time, the kinatic energy spectrum changes quickly 
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and r«-aatabliahas tha isl noda at tha ona with tha largaat 

atiount of kinatic anargy. By t«.152«tha aolution lattlaa down 

from ita rapid oacillations. Tha pariod of tha anargy apac- 

trum ia t e.l528, however, tha ptriod of tamparatura and valo- 
P % % 

city ia 2t^. Wa have fourd that T(f*-t_)»-T(t) and v(ttt^)* 

P p - p 

-v(t^t-). If we aaauaa that tha charactarlatie velocity of 

P in 

tha fluid ia C2 K£ (i*l, ,than wa can aatioata 

tha addy turnover time, t^, to ha C2 KE (£«!, r*5) 
or 0.022. Tha period of tha apactrum, t^ ia 6.95 t^. Wa 
have repeated tha calculation with with tha 

viacoaity of tha i*3 oodaa (but not tha t«l or 2 modaa) incraaaad 
by 101. With tha enhanced viacoaity tha solution ia steady- 
atata. Whan wa increased tha thermal diffusivity of tha I«3 
modes by lOt, tha solution remained periodic in time. 

With ^outoff*^* solution is both time dependant and 
aperiodic. Tha solution wanders between a normal 
and an inverted state. In tha normal state KECisD^lO^ and 
KE(4s2>2lO. In the inverted state KZ(i»l):iC^ and KE(A82)tlO^. 
The time dependence of the solution ia raminlacant of tha ]r.:kr.nar 
in which a Lorenz solution wanders between two strange attrac- 
tors. We have not attempted to determine whether there are 
fixed points in the ^^utoff*^ equations of motion. 
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IV. DISCUSSION 

W« h«v« oalculAttd solutions to tht squstions of convsc- 
tion by sxp«nding tht horizontal structurs in s ssriss of aodst. 
Choosing ths numbsr of aodss to bs rstsintd in tht solution is 
squivalsnt to fixing tht spatial rtsolution in tht horizontal 
dirtction. Kttping tht numbtr of aodts fixtd, wt found that as 
tht Bayltigh numbtr incrtasts, tht solution changts in tiat from 
sttady-statt to ptriodic, and than to aptriodic. Alttrnativtly, 
wt havt found that by Vttping tht Rayltigh numbtr fixtd 
at 10^ and dtcrtasing tht numbtr of modts, tht solution changts 
from sttady-statt t to ptriodic » to aptriodic. In tht txtrtmt 
cast whtrt tht txpansion is limittd to modts of only ont hori- 
zontal wavtnumbtr» tht solution gots to a sttady-statt singlt 
modt. From our obttrvations of tht timt-indtptndtnt solutions 

4 

with Rj*10 , wt havt found that truncating tht horizontal 
txpansion rtsults in: 1) alttring tht kinttic and thtr- 
mal spifctra by incrtasing tht amplitudts of tht high wavtnuabtr 
modes t 2) maxing tht mtan ttmptroturt graditnt mort isothtrmal 
and thtrtby lowtring tht ctntral ttmptraturt, and 3) dtcrtasing 
tht ratt at which tht ttmptraturt varianct is produetd in tht fluid. 
Wt havt shown that if tht truncation is too stvtrt^tht thermal 
varianct spectrum will btcomt inverted, with tht high wavenumber 
dissipation modts having mort energy than tht low wavenumber 


production modes. Tht thermal varianct inversion dots not de- 
stroy tht time- independent property of the fluid. We 
have shewn in one example that if the truncation is severe enough 
that the kinetic energy spectrum is inverted, then the solution 
becomes time dependent. T.his suggests that we should be cautious 
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in numerically computing points of bifurcations in the flow 
as a function of Rayleigh number alone. The points of bifur- 
cations become curves of bifurcations when plotted both as a 
function of Rayleigh number ^uld limit of spatial resolution. 
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FIGURE CAPTIONS 
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Figure 1 - The kinetic energy for the £=1 mode as a function of 
radius calculated with (solid line) and 

(broken line). The higher kinetic energy in the single- 
mode calculation allows more kinetic energy to be viscously 
dissipated and compensates for the inability of the single- 
mode calculation to lose energy by cascading. 

Figure 2 - Same as figure 3 with the temperature variance of 
the £=1 mode plotted as a function of radius. 


Figure 3 - Q = KE(£ = 2, r=0. 5)/KE(3.r=0. 5) as a function of 
^cutoff* Truncation causes the high wavenumber modes of 
the kinetic energy spectrum to become anomalously large. 

By extrapolation, it appears that when ^g^^off'^* 

meaning tliat the kinetic energy spectrum has become inverted. 

Figure 4 - The kinetic energy calculated with 

r=0.5 for the £=1, 2, and 3 modes as a periodic function 
of time. At t=.0603 kinetic energy inverts so that 
KE(£=1, r=.5) < KE(£=3, r=.5). 
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